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Abstract 

In this paper we report calculations of some pseudospin textures for bilayer 
quantum hall systems with filling factor = 1. The textures we study are 
isolated single meron solutions. Meron solutions have already been studied 
at great length by others by minimising the microcopic Hamiltonian between 
microscopic trial wavef unctions. Our approach is somewhat different. We 
calculate them by numerically solving the nonlinear integro -differential equa- 
tions arising from extremisation of the effective action for pseudospin textures. 
Our results can be viewed as augmenting earlier results and providing a basis 
for comparison. Our differential equation approach also allows us to dilineate 
the impact of different physical effects like the pseudospin stiffness and the 
capacitance energy on the meron solution. 

I. INTRODUCTION 

Recent years have witnessed the demonstration of topologically nontrivial spin textures 
like Skyrmions as physically important excitations in quantum Hall systems. It has been 
argued in a body of very interesting papers that such Skyrmions are not mere exotic theo- 
retical curiosities but will actually be the lowest energy excitations in some situations ( such 



as for example filling factor u = 1 and small Zeeman coupling) [Q. Furthermore, experimen- 
tal evidence indicating the presence of such Skyrmionic excitations has also subsequently 
emerged 0. Skyrmions of topological charge more than unity have also been studied 0]. 

These results were basically for mono-layers of electron gas in strong magnetic fields, 
but with spin degrees of freedom of the electrons treated dynamically. Meanwhile, the 
remarkable discoveries of the quantum Hall effect, originally found in single two dimensional 
electron layers, have also been extended to bilayer systems , thanks to the development 
of techniques for growing GaAs heterostructures containing two separated layers of two- 
dimensional electron gas (see for example references @]) . On the theoretical front, a large 
body of work has already been done on bilayer systems. An extensive list of references to 
this literature has been given in the lucid review of this subject by Girvin and MacDonald 
and in the paper by Moon et al 0. 

To simplify matters, one may begin analysing bilayer systems by assuming , to start 
with, that the electrons are spinless, i.e. that the spins of the electrons are fully polarised 
by the strong magnetic field and frozen as degrees of freedom. The idea would be to later 
incorporate the complications brought about by spin degrees of freedom after the spinless 
problem is understood. However, even when real physical spin is suppressed, it is helpful 
to view electrons in a bilayer system as carrying a "pseudospin" in addition to their spatial 
coordinates on the plane. The notion of such a pseudospin arises from mapping the bilayer 
spinless problem into a monolayer problem with spin 0. In this mapping the layer degree 
of freedom is specified by a two-component pseudospinor whose up and down components 
refer to the amplitude for the electron to be in the first and second layers respectively. 
Such a mapping allows one to borrow for bilayer systems, the rich body of insights and 
results available from single layer systems with real spin. For instance, just as in case of 
electrons with real spin the combination of the Pauli principle and Coulomb interactions can 
lead to a ferromagnetic ground state, one can similarly expect the pseudospin of bilinear 
systems to also be fully "polarised" under corresponding circumstances. Such psuedospin 
magnetisation amounts to spontaneous phase coherence between the two layers, with many 
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interesting consequences , discussed in the literature 0,0. 

Going beyond the ground state, the relevance of Skyrmions to systems with real spin 
has in turn prompted studies of similar topological excitations in spinless bilayer systems, 
but now involving pseudospin (See references [Qand references given therein). Because 
of interplane-intraplane anisotropy in pseudospin stiffness in bilayer systems, as well the 
capacitance energy of carrying unequal charge density in the two layers, the effective Action 
governing pseudospin enjoys only U(l) symmetry of rotations about the z-axis (the direction 
perpendicular to the x-y plane of the layers). Finiteness of energy requires that asymptoti- 
cally the pseudospin must lie on the easy (x-y) plane. The basic topological excitations in 
that case are the so-called merons which are vortices in pseudospin with a winding num- 
ber of one-half, and meron anti-meron bound pairs. Such a pair is topologically equivalent 
to Skyrmions and carries unit winding number. (For an introduction to such topological 
excitations, their winding numbers, requirements due to finite energy, etc. see reference 0.) 

In this paper we obtain single meron solutions in bilayer systems with filling fraction 
1/ = 1, for a range of different values of system parameters (such as the inter layer separation 
d) by numerically solving the integro differential equation obeyed by them. Even though a 
single meron by itself has a logarithmically divergent energy in the infinite volume limit it is 
nevertheless an object of interest. Physical meron anti-meron solutions which do have finite 
energy, are often treated as a bound pair of two extended interacting objects, namely the 
individual merons each with a " size" and a meron core energy . Merons are not rigid bodies 
but rather solutions of nonlinear wave equations; the notion of giving each of them a fixed 
size, shape and energy, although only approximate when they appear as a part of a larger 
solution containing other merons, is nevertheless useful in such Soliton physics. Single meron 
solutions allow us to obtain their size, shape and individual energy in an unambiguous way. 

Single meron and bi-meron solutions have already been extensively studied in a body of 
papers by Girvin, MacDonald and co-workers |Q ||10|| and . These calculations are based 
on optimising microscopic wavefunctions with respect to the microscopic interaction Hamil- 
tonian. The emphasis is on determining the optimised energy of bimeron configurations and 
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extracting individual meron core sizes and energies from them. Our work here rehes heavily 
on the advances already made in these papers, and is to be viewed as something which will 
hopefully augment their calculations. We will obtain explicit meron solutions through a 
somewhat different route , by numerically integrating the nonlinear differential equations 
derived from the effective action obeyed by the pseudospin magnetic moments. Meron pro- 
files are of some interest in their own right. [It is for instance not widely realised that unlike 
Skyrmions which exist even in the pure nonlinear sigma model (NLSM) limit, that model 
will not permit exact meron solutions. Meron solutions require departure from the NLSM 
model, which is provided by the capacitance energy and other anisotropies. We will explain 
this point further later on.] We will obtain our meron profiles for a range of values of the 
layer distance and study the way in which the solution is affected by the addition of leading 
corrections to the NLSM effective action. We also compare our results for the meron size 
and core energy with those obtained from other methods. 

II. RECAPITULATION OF THEORETICAL PRELIMINARIES 

The integro-differential equation which we will integrate to get our meron solutions is 
obtained by extremising an effective action which has already been derived by Moon et al 
starting from the basic microscopic physics. See also Ezawa |Tl|]. For completeness , and 
to set up the notation let us quickly recall their result . The pseudospin texture of a state is 
described by a classical unit vector m(r) which gives the local direction of the pseudospin. 
Here r is the coordinate on the x-y plane carrying the layers, while the magnetic field B 
is along the z-direction. The fully polarised "ferromagnetic" ground state |\E'o) corresponds 
to m pointing everywhere in the same direction, say, along the x-axis. Using this as the 
reference state, any other state with some arbitrary texture m(r) is given by performing a 
local pseudospin rotation on this uniform ground state : 




(2.1) 



Here the pseudospin rotation operator is, for infinitesimal rotations. 
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O = J d^x {xx m(f)) . S{r) (2.2) 

where S is the local pseudospin operator given hj S = ^ipl^ ipp and ipa is the second 
quantised electron field operator with a denoting the pseudospin (layer) index. The effective 
action I[m\ for any given pseudospin texture m(r) has already been derived by Moon et al 



0] starting from the basic microscopic physics. See also Ezawa [|rT|. The microscopic 
Hamiltonian is 

H = Hk + Hint (2.3) 



where 



and, 



Hk = ^jzjd'x (2.4) 
Hint = \\ d^xd^x' 5p„(x) K/3(x - x') %(x') (2.5) 



and 6pa is the local density operator minus the mean density in the ath layer. The interaction 
potential written here as a 2 x 2 matrix in pseudospin (layer) space is due to Coulomb forces 
between electrons in the same or different layers and is given by 

y„;3(x-x') = , (2.6) 
e^(x-x')2 + d2(l-M 

where dis the layer separation. 

The effective action I[m\ is just the expectation value of the microscopic Hamiltonian 
in the state |m(r*)), but with all operators projected on to the lowest landau level (LLL). 
The detailed calculation is given in ref P| and leads to a nonlocal action functional of m(r). 
Upon doing a gradient expansion , valid for slowly varying (long wavelength) textures, the 
leading terms are 

/(m) = J d^r [lp^(Vm,)=^ + ip^((Vm,)2 + {VmyY) + (3ml^ + C[m] (2.7) 

where 



and terms of order and higher have been dropped. The constants appearing in this 
effective action are given in terms of microscopic interaction parameters. The constants pA 
and pe are pseudospin stiffness parameters whose physical origin is the exculsion principle 
(Hund's rule) mentioned earlier. They are given by 

PA = (t^) / dkk^ exp{^) 

^ ^32^) i exp(^) (2.9) 

where = 27re^/(eA;) and = ea;p(—/cd)(27re^/eA;) are the Fourier transforms of the 
Coulomb interactions between electrons in the same and different layers respectively. All 
distances (and inverse wave vectors) are in units of the magnetic length /. The (3ml term 
represents the so-called capacitance or charging energy needed to maintain unequal amounts 
of charge density in the two layers. Recall that the z-component of pseudospin represents 
the difference between the densities in the two layers. The constant (3 is given by 

wherey,^ = |(y/-lf^). 

III. THE FIELD EQUATION AND ITS SOLUTIONS 

The nonlinear partial differential equations that have to be obeyed by all slowly varying 
spin textures, whether they be single merons or multi merons is obtained by extremising this 
effective action I[m] with respect to m{r). In doing this extremisation one has to respect 
the condition that m is a unit vector. We handle this by writing m in terms of independent 
fields rUz and a such that 

m. = [Jl — ml cosa, Jl — mlsina, mz) (3.1) 
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As one can see a = tan~^{^) is the azimuthal angle of the vector m in the x-y plane. 
Upon rewriting J[m] in terms of the independent fields m^(r) and a;(r*), and minimising 
with respect to them we get the field equations 

(3.2) 

and 

V.[(l -m^)Va] = (3.3) 

We will numerically solve these coupled equations for the variables and a . In the 
case of the Skyrmion and bi-meron solutions people have found it helpful to begin with the 
Isotropic Nonlinear Sigma model (NLSM) limit of the effective action and field equations, 
whose analytical solutions are exactly known. These exact solutions offer a starting point 
for setting up trial functions for numerically solving the full equations. In the case of the 
meron however, there is no exact single meron solution to the NLSM. To see this let us 
briefly consider the NLSM limit. 

A. Absence of single-merons in the NLSM limit 

The SU(2) invariant NLSM limit is obtained when the layer separation vanishes {d/l — >■ 
0). At the level of the microcopic Hamiltonian one can see that it is rotationally (SU(2) ) 
invariant in pseudospin in the limit of zero separation. The potential energy in ( |2.5| ) and 
( p.6| ) is layer independent when d = 0, and in the LLL approximation the kinetic term 
( p.4| ) involving the magnetic fleld is just a constant . This symmetry is also reflected in 
the effective action I[m\. With(i= 0, we see from their deflnitions that pA = Pe while 
/? = = C[m]. Then the effective action in ( p^ ) reduces to just the NLSM. The equation 
for TJiz in turn reduces to 



while |3.3| remains the same. The NLSM has been fully solved and its solutions are most 
conveniently described in terms of the complex field w(z) defined by 

where z = x+iy. Any analytic function w(z) will be a solution of the NLSM. In particular, 
the single Skyrmion solution is w{z) = |, which corresponds in polar coordinates (r, 6) to 
the circularly symmetrical function with 

r2 - A2 

iT^z = o,nd a = 9 (3.6) 

_l_ \^ 



It can be checked that (p.6|) satisfies equation ( p.4|) 

Now , isolated meron or anti-meron solutions may be defined by the boundary conditions 

■mz{r = 0) = ±1 

and rriz = ; rrix + iniy = exp{±i9), as r oo (3.7) 

Solutions with such boundary conditions may or may not exist for a given equation. In fact 
we will argue that the basic NLSM has no meron solutions. But the NLSM does have a 
Skyrmion solution of the form (|3.6|). The notion of a meron stems from the portion of the 



Skyrmion function (|3.6|) from r = to r = A. In this portion the circle of radius A in the z- 
plane is mapped onto the unit circle in the w-plane. The spin begins at r = with = — 1 
starts tilting towards the x-y plane, and by r = A is lying entirely on the x-y plane, pointing 
radially outward. The topological number density when integrated over this portion, is one- 
half. Sometimes this portion is called a meron, i.e. half a Skyrmion. However, it must be 
emphasized that this object corresponds only to a piece of the full Skyrmion function and 
is not a solution to the field equations in the full x-y plane. In contrast an "isolated single 
meron solution" to the field equation, while qualitatively similar, should span the whole x-y 
plane and not just a portion of it. That is, it should satisfy the boundary condition ( p.7| ) 
at infinity and not at some finite r = A. Such a solution does not exist. The easiest way 
to understand this is to note that there exists no analytical function w(z) which maps the 
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z-plane into a unit circle in the w-plane. Equivalently, if one tried to solve the differential 
equation |3.4j starting with = —1 at r=0 and any positive slope, one will always end up 
with a Skyrmion, i.e. with = +1 at r = oo. Note that this is not directly related to 
the fact that an isolated meron ,if it existed, would have a logarithmically divergent energy. 
Solutions can exist to diferential equations with some boundary conditions, even if they have 
infinite value for the energy functional. But in the case of merons, no such solution exists 
for the NLSM. 



B. Meron solutions when d 7^ 

When the layer separation d is non-zero, then the full effective action (|2.7| ) comes into 
play and we have no longer just the NLSM. The arguments of the preceding subsection do 



not hold and now meron solutions can exist. See Affleck for arguments show that in 
fact they will exist when (3 in eq ( |3.2| ) is non-zero. 

We will look for spherically symmetric solutions with a(r, 0) = 9 which will satisfy 
eq(|3.3|) provided depends only on the radius r . Then eq (p. 21) requires that m^(r) satisfy 
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2(jm-^ - pA[Orm^ + -drm^) - pe\ —. ^ H ^ h 

r \ (1 — m'^^Y 1 — mi, J 

The parameters pA, Pe and /5 in the equation are calculated for each value of the layer 
separation d from their definitions ( |2.9| and p.lO| ) . 

As it stands ( p.8| ) is a nonlinear integro-differential equation. It is the nonlocal integral 
in the last term which represents the most difficult part of this equation. To start with let 
us neglect this term . We are then left with a nonlinear ordinary second order differential 
equation given by 

2/?m, - Pa(<9X + -drm,) - pe , ^ Z + \ I + ^ (3-9) 



r V (1 — mlY 1 — ml / 

We will first solve this equation and then later on return to the full equation ( p.8|) 
including the last integral term. To solve (|3.9| ) numerically we merely start with a boundary 
value and slope for at r = , obtain d^mz using the equation ( |3.8| ), use this to get the 
value and slope at the next point in r , and thus proceed towards r — ^ oo. One can see that 
the aymptotic behaviour of mz{r) must be 

mz = 1 — ar^ + O(r^) as r 
and mz ~ exp{—kr) as r ^ oo (3.10) 

where k = and a is a constant for a given layer separation d. The absence of an 

0(r) term in m^(r ^ 0) follows from the requirement that our solution not have a cusp at 
r = 0. One may be concerned that the last term (proportional to pe) in eq (|3.9|) appears 
to have pieces with an singularity at r = 0. It must be remembered that (1 — m^)~^ 
also behaves as r~^. However it can be checked that the requirement ( p.lO|) ensures that 



the singularities in the different pieces cancel one another, regardless of the value of the 
constant a . This constant a is adjusted so that exponentially falls to zero as r — oo as 
required. 
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IV. RESULTS AND DISCUSSION 



The resulting solutions of eq (|3.9| ) for mz{r) are plotted in figure 1 for a range of 6 
different values oi d = 0.24 , 0.50 , 0.625 , 0.78 , 1.00 , 1.20 and 1.50 respectively in units 
of the magnetic length / . Each solution falls monotonically to zero as r increases with an 
exponential fall as r ^ oo as required by (|3.10| ). Comparing solutions for different values of 
the interlayer distance ci, we see that the width (mean radius) of the solution monotonically 
decreases with ci, the outermost curve shown in fig.l corresponding io d = 0.24 and the inner 
most to d = 1.50. The physics behind this is the minimisation of the capacitance energy 
f3ml , which would tend to bring down to zero as quickly as possible, in contrast to the 
gradient (stiffness) terms would tend to reduce as slowly as possible from its starting 
value of unity at r = 0. As d increases so does /3 , while pe falls with d and pa is independent 
of d (see equations ( |2.10|) and p.lO|) . Therefore the relative importance of the capacitance 
energy to the stiffness terms increases with d, and hence the solutions fall to zero faster. 

In fig. 2 we plot the meron core radius Rmc, which we define following the literature |10 



as the radius at which equals 0.1. As noted above this radius falls monotonically with 
d. In fig. 2 we show our values of Rmc for 8 different values of d, while for comparison the 



results of Yang and MacDonald (extracted from fig. 2 of ref. |jTO[|) are shown in the form of 
a curve. They consider their results to be reliable in the region 0.5 <d< 1.2. In this range 
our values of Rmc are of the same order of magnitude as theirs but about 30 to 40 percent 
smaller . The qualitative behavior as a function of d is very similar. 

The energy of the meron, i.e. the value of the effective action for the meron solution 
will diverge logarithmically if evaluated over the entire plane because of its nonvanishing 
azimuthal gradient. However the notion of a "meron core energy" applies to the energy of 
the meron function evaluated upto r = Rmc- In fig-3 we plot Emc as a function of d. 

Since single merons have , strictly speaking, divergent energy, physical finite energy 
excitations involving merons must come in the form of meron-antimeron pairs, which do 
have finite energy . Using Emc one can make a phenomenological estimate of the energy 
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Emp of such meron pairs as pointed out in references and [0. Following them we can 
approximately write the energy of a pair separated by a distance R in the form 

Emp = 2Emc + 2npEln{R/R^c) + -r-B (4-1) 
The optimal value R* of the meron separation can be obtained by minimising this energy 

2 

with respect to R to get R* = j^:^^^- Substituting this value of R* into the above equation 
we get 

Emp = 2E^, + 2npE(l + H—^-—)] (4.2) 

Table 1 displays the values of meron core radius, core energy, the merom pair separation 
and energy as obtained by us, for different values of d. These may be compared to values 
obtained by others using microscopic Hartree Fock methods. We see that for the lower 
values od d , such as 0.24 and 0.5, the meron separation R* is either smaller or of the same 
order as the core radius Rmc- The approximation of two well separated merons is clearly not 



good at these values of d. Equations such as ( [4.1[ ) should be used only at layer separation 
d beyond 0.6. This is in agreement with the findings of Yang and MacDonald |10[ . 

Finally we include the non-local integral (the last term) in the original field equation 



^ which we had neglected so far. It arose in this equation as the functional derivative 



^^^^ of the term C[m\ in the Effective action (|2.7| ). The presence of this term not only 



Sniz 



makes the equation more difficult to solve , but also alters to some extent the nature of the 
solution. Notice that this term behaves for large r as ^ where C oc / d'^r'V'^mz{r'). This ^ 
behavior will not permit an exponential fall off of as had happened in the solutions in 
the absence of this term , shown in figure 1. Instead the solution can be expected to fall to 
zero asymptotically only as some inverse power of r. Further, the value of the constant C 
for most values of d of interest is positive. For large r the solution for m^., may therefore 
approach zero from below. 

These expectations are supported by the numerical solution of the integro-differential 
equation (|3.8| ). We solve it by the following iterative procedure. Start with the solution m° 
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of eq( p.9|) , where the above-mentioned integral has been altogether neglected. Evaluate the 
integral by inserting in its integrand. Multiply the resulting integral by a prefactor a = 
0.1 , add it to eq.(|3.9|), and solve the equation again to obtain m°. Solve the equation again 
after inserting this into the integral, call the solution m^°, and repeat this procedure ten 
times until you get m^,°". Since lOa = 1 this should yield a good approximation to the full 
equation ( ^78|) . One could have used this iterative procedure in smaller steps , but we find 
the step size of a = 0.1 is small enough for the procedure to converge smoothly . 

In figure 4 we compare the solution so obtained for the full equation ( p.8|) with the 
"unperturbed" solution (of the simpler eq ( p.9|) discussed earlier), for two values of d equal 
to 0.625 and 1.20. One can see that the full solution reaches zero, and then after turning 
negative, returns to zero asymptotically but more slowly than the unperturbed solution. In 
other words the pseudospin, which starts out pointing along the z-axis at r=0, gradually falls 
towards the x-y plane, lies entirely on the x-y plane at some finite radius, then starts bending 
towards the negative z-axis , then reverses its behavior to finally asymptotically come back 
to the x-y plane at infinity. We see no objection to such a meron solution. But the notion 
of a core may be less meaningful now. One could formally evaluate the "meron core radius 
RmC " of this new solution defined once again as the radius at which = 0.1, and a 
corresponding core energy. But with the solution continuing to be substantially nonzero at 
radii quite a bit larger than such an RmC , the concept of a core is less significant for such 
slowly falling solutions. 
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Table 1: Meron core radius (Rmc), meron-antimron separation (i?*),meroii core energy ( Emc 
) and meron pair energy (Emp) as a function of layer separation. The unit of energy is ^ 
and the unit of length is / 



d 


Rmc 


R* 


E 


Emp 


.2400 


5.07276 


2.3133 


0.145487 


0.314187 


.5000 


3.07692 


3.52110 


0.105005 


0.290580 


.6250 


2.70270 


3.9868 


0.101588 


0.290250 


.7810 


1.99584 


4.88985 


0.084911 


0.266760 


1.000 


1.41372 


6.42600 


0.067163 


0.232140 


1.200 


1.14761 


8.15756 


0.057474 


0.205750 


1.500 


0.873181 


11.3682 


0.042215 


0.163039 
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FIGURES 

FIG. 1. The solution mz{r) of equation ( |3.9D . The six curves correspond to layer separation d 
equal to 0.24 , 0.50 , 0.625 , 0.78 , 1.00 , 1.20 and 1.50 respectively in units of the magnetic length 
/. The width of the curves decreases monotonically as as d increases, with the outermost curve for 
d = 0.24 and the innermost to d = 1.50 

FIG. 2. Meron core radius Rmc for 8 different values of the layer separation d , all in units of 
/. The solid curve gives the results obtained in ref.lO 

FIG. 3. Meron core energy Erac a function of layer separation 

FIG. 4. mz{r) versus r. The continuous lines give the solution of eq( |3.8f ) which includes the 
integral term and, for comparison, the dotted curves depict the solution of ( |3.9D which does not 
include this term. Figure (a) corresponds to d = 0.625 and figure (b) to d = 1.2 
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